No-Core MCSM calculation for 10 Be and 12 Be low-lying spectra 
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The low-lying excited states of 10 Be and 12 Be are investigated within a no-core Monte Carlo Shell 
Model (MCSM) framework employing a realistic potential obtained via the Unitary Correlation 
Operator Method. The excitation energies of the 2f and 2 J states and the B(E2; 2^ 2 — > 0^ s ) 
for 10 Be in the MCSM with a standard treatment of spurious center-of-mass motion show good 
agreement with experimental data. Some properties of low-lying states of 10 Be are studied in terms 
of quadrupole moments, E2 transitions and single-particle occupation numbers. The E2 transition 
probability of 10 C, the mirror nucleus of 10 Be, is also presented with a good agreement to experiment. 
The triaxial deformation of 10 Be and 10 C is discussed in terms of the B(E2) values. The removal 
of the spurious center-of-mass motion affects differently on various states: for instance, negligible 
effects on the 2f and 2% levels of 10 Be, while significant and favorable shift for the 1]~ level. It 
is suggested that the description of 12 Be needs a larger model space as well as some other higher 
excited states of 10 Be, as an indicator that these are dominated by intruder configurations. 



PACS numbers: 21.60.Cs,21.60.De,21.60.Ka,27.20.+n 



I. INTRODUCTION 

One of the major goals in nuclear physics is to un- 
derstand the structure and reactions of nuclei starting 
from realistic nuclear interactions. Besides the challenge 
of solving the nuclear many-body problem, this endeavor 
is complicated by the fact that our understanding of the 
nuclear force is not complete yet. At present, there are 
two ways to construct an accurate representation of nu- 
clear force. One can construct a two-body potential phc- 
nomcnologically by fitting experimental data on nucleon- 
nucleon (NN) scattering, as it is done in the Argonne 
V18 potential Q, the CD-Bonn potential and the Ni- 
jmegen potentials [||. Alternatively, consistent two- and 
many-body interactions can be constructed in the frame- 
work of chiral effective field theory using the symmetries 
and the effective degrees of freedom of low-energy QCD 
as a guiding principle. The chiral N 3 LO potential is such 
an accurate charge-dependent nucleon-nucleon potential 
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constructed at fourth order of chiral perturbation the- 
ory By using these realistic nuclear interactions, 
ab initio nuclear many-body calculations have been per- 
formed in the last decade. In Green's Function Monte 
Carlo (GFMC) calculations the exact ground-state wave 
function is calculated by treating the many-body Green's 
functions in a Monte Carlo approach The GFMC 
calculations of light nuclei up to 12 C with the Argonne 
interaction reproduce the experimental nuclear binding 
energies and radii as well as the spectra. Another ab 
initio approach for nuclei up to A=14 is the No-Core 
Shell Model (NCSM) All nuclcons arc treated in 
a large number of shell-model basis, providing similarly 
successful description of light nuclei. 

However, the straightforward application of those re- 
alistic interactions in nuclear many-body calculations 
is still difficult due to the strong short-range repulsion 
which generates strong correlations in the nuclear many- 
body state. The Unitary Correlation Operator Method 
(UCOM) is one of the methods to tackle this problem by 
introducing a unitary transformation such that the trans- 
formed many-body states contain the information on the 
dominant correlations in nuclear many-body system []~3l — 
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Il5| . In the UCOM approach two unitary transformation 
operators are defined, a central correlation operator and 
a tensor correlation operator, which correspond to two 
most important correlations: the central correlations in- 
duced by the strong short-range repulsion and the tensor 
correlations, respectively. Through a unitary transfor- 
mation of the Hamiltonian, a soft phase-shift equivalent 
two-nuclcon interaction can be obtained. This UCOM 
potential can be used in various kinds of many-body cal- 
culations, such as no-core shell model calculations flU - 

Q. 

In the shell model calculations, the direct diagonaliza- 
tion of the Hamiltonian matrix in the full valence- nucleon 
Hilbcrt space is difficult, as the dimension of such a space 
becomes larger and larger when one moves from light nu- 
clei to heavier nuclei. Much effort of truncation frame- 
works to full shell-model calculation has been directed, 
e.g. in Refs pjl H(J As another way to overcome 
this difficulty, the stochastic approaches have been in- 
troduced. Among them, the Shell Model Monte Carlo 
(SMMC) method has been successfully proposed f22j . 
Nevertheless, the SMMC is basically suitable for the 
ground state and thermal properties, and suffers from 
the so-called "sign problem". As a completely differ- 
ent approach, the Quantum Monte Carlo Diagonalization 
(QMCD) method has been proposed for solving quantum 
many-body systems with a two-body interaction [23U2(| . 
The QMCD can describe not only the ground state but 
also excited states, including their energies, wave func- 
tions and hence transition matrix elements. Thus, on the 
basis of the QMCD method, the Monte Carlo Shell Model 
(MCSM) has been introduced [13] for nuclear shell model 
calculations (2^-[13- An extrapolation method in the 
Monte Carlo Shell Model has been proposed very re- 
cently [33[. The applicability of the MCSM to a sys- 
tem beyond the current limit of exact diagonalization is 
shown for the pf + gg /2-shell calculation by assuming a 
core in their work. It is then of a certain importance and 
interest to apply the MCSM to ab initio calculations of 
light nuclei. As the MCSM has never been used in ab 
initio calculations, we start with straightforward calcu- 
lations by taking conventional MCSM method and code 
which have been used for many shell-model calculations 
for medium-mass nuclei. We shall present, in this paper, 
how such ab initio calculations work. We note that the 
MCSM method is being revised in parallel, and outlines 
of such revisions and future directions can be found in 
Refs. [HI, [H[ . The results to be shown in this paper will 
play a key role in judging as to whether one should move 
ahead to more systematic calculations with the revised 
method or not. 

The MCSM calculation is performed without a core 
to make it ab initio. In Section 2, we will introduce the 
theoretical framework of the MCSM and explain the gen- 
eral procedure of the Monte Carlo Shell Model method. 
4 He, which is investigated in the framework of the shell 
model and the MCSM by using the UCOM potential, is 
discussed in Section 3 as the numerical check. Study of 



structure and low lying spectra for Beryllium isotopes 
appears in Section 4. In Section 5, the conclusion with a 
summary and description of future direction for research 
in this field is given. 



II. MONTE CARLO SHELL MODEL 
CALCULATION 

The main idea of the MCSM is to diagonalize the 
Hamiltonian in a subspacc spanned by the MCSM ba- 
sis states, which are generated in a stochastic way. 

We begin with the imaginary-time evolution operator 

e-" H , (1) 

where H is a given general (time-independent) Hamilto- 
nian and (3 oc T^ 1 is a real number with T being anal- 
ogous to a temperature. If this operator in Eq. (|T|) acts 
on a state |\p( )), one obtains 



-/8#|$(°)) 



-/3E, 
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(2) 



where Ei is the z-th eigenvalue of H, is the corre- 
sponding cigenstate and c, its amplitude in the initial 
state: 



(3) 



For (3 large enough, only the ground and low-lying states 
survive. But the actual handling is very complicated for 
H containing a two-body (or many-body) interaction. 

The Hubbard-Stratonovich (HS) transformation [HI, 
[37| can be used to ease the difficulty mentioned above. 
We then move to the formula 



|*(<r)> oc e -^M|*(°)), 



(4) 



where h(a) is a one-body Hamiltonian obtained through 
the HS-transformation and a is a set of random numbers 
(auxiliary fields). The right-hand-side of this relation can 
be interpreted as a means to generate all basis vectors 
needed for describing the ground state and the low-lying 
states. For different values of the random variable, <r, one 
obtains different state vectors, |$(tr)), by Eq. (J4j> . These 
vectors are labeled as candidate states and selected as 
MCSM basis by a procedure of energy comparison (2?J . 

During the MCSM generation of the basis vectors, 
symmetries, e.g. rotational and parity symmetry, are 
restored before the diagonalization as more basis vectors 
arc included. All MCSM basis states are projected onto 
good parity and angular momentum quantum numbers 
by acting with the corresponding projection operators. 
We diagonalize the Hamiltonian in a subspace spanned 
by those projected basis vectors. The number of the 
MCSM basis states is referred to as the MCSM dimen- 
sion. The basis generation process for general cases is 
outlined in Ref. [27| . 

As more than one major shell is included in the MCSM 
calculation, the spurious center-of-mass motion must be 
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accounted for. The Lawson's prescription is adopted to 
suppress the spurious center-of-mass motion in good ap- 
proximation for major shell truncation [38[. The total 
Hamiltonian then consists of intrinsic and center-of-mass 
parts as, 



H = Hi, 



where H int . 
defined by 

He. 



is the intrinsic Hamiltonian. The H c 



-MAu 2 R 2 - -hu, 
2 2 



(5) 



is 



(6) 



2AM 

where R and P are the coordinate and momentum of 
the center of mass, respectively. In general, by taking 
sufficiently large values of /? c .m., spurious components are 
suppressed for the low- lying eigenstates of H' . 

III. RESULTS FOR 4 He 
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FIG. 1. (Color online) The ground-state energy of 4 He as 
a function of harmonic oscillator frequency hw in different 
model spaces (e m ax from 2 to 5) calculated in NuShell us- 
ing the Vucom(N 3 LO) potential. The symbols correspond to 
e ma x from 2 to 5. The straight line is experimental value. 



In this section we discuss the interactions and model 
spaces used for the no-core MCSM and provide some 
benchmark calculations for the 4 He ground state. The 
model space of the MCSM is spanned by a harmonic os- 
cillator basis truncated with respect to the unperturbed 
single-particle energies e max = 2n + I. We use UCOM- 
transformed realistic two-nucleon interactions as input 
potential. In addition to the standard UCOM interac- 
tion derived from the Argonne V18 potential, which has 
been used in a series of applications in various many-body 
methods we adopt a new UCOM potential based on 
the chiral N 3 LO two-nucleon interaction of Entem and 
Machleidt 0,0. These UCOM potentials are labeled as 
Uucom(AV18) and Vu C om(N 3 LO), respectively. In both 
cases the UCOM correlation functions are determined 
through an energy minimization in the two-nucleon sys- 
tem with a constraint on the range of the even spin-triplet 
tensor correlator [l5| . We neglect Coulomb interaction in 
all of our calculations throughout this work for simplicity. 

As an example of the UCOM potential, we perform 
a straightforward shell model diagonalization within a 
harmonic oscillator basis without a core. In this shell 
model calculation, we employ the NuShell code developed 
by B.A. Brown et al. [39| and use the Vucom(N 3 LO) 
potential. 

The ground-state energy for 4 He as a function of the 
harmonic oscillator frequency hu in various model spaces 
characterized by the oscillator basis cut off parameter 
e m axj is shown in Fig. Q] The shell- model description of 
many-body correlations depends in general on the size 
of the model space. The experimental value is shown as 
a black line. The ground-state energy for small model 
spaces, e.g., e max =2, shows a sizable dependence on hui. 
By increasing the size of the model space, the ground- 
state energy is lowered and dependence on hui is re- 
duced. The ground-state energy varies by about 1 MeV 
for a range of oscillator frequencies hui from 24 MeV to 
52 MeV. There is still about I MeV difference between 
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FIG. 2. (Color online) The ground-state energy of 4 He as 
a function of the harmonic oscillator parameter hui using 
UCOM potential Vucom(AV18) in e max =l (black), 2 (blue), 
and 3 (red) model space. The open symbols indicate results 
by the conventional direct diagonalization method in the m- 
scheme. The closed symbols indicate the MCSM results. 



the e max =5 result at hui=32 MeV and the experimen- 
tal ground-state energy. Evidently, still larger values of 
emax are needed to reproduce the experimental binding 
energy of 4 Hc, which the Vucom(N 3 LO) interaction is 
approximately adjusted to. At present, the conventional 
shell model calculation is performed only up to e max = 5 
model space due to the limitation of Nushell code. How- 
ever, the convergence is significantly better than for the 
bare realistic nuclcon-nuclcon potential. The compari- 
son between bare interaction and the UCOM interaction 
in the shell model calculation for 4 He can be found in 
Ref. HU. 

We now compare the full shell-model and the MCSM 
results for the ground-state energy of 4 He using the 
Hjcom(AVI8) potential. Figure [2] shows the ground- 
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FIG. 3. The center-of-mass motion energy and the ground- 
state energy (inset) of 4 He as a function of Lawson's prescrip- 
tion parameter /3cm. defined in eq. ((5}. The model spaces are 
e ma x=l, 2 and 3. The oscillator parameter hw is 28 MeV. The 
potential Vtjcom(N 3 LO) is used. 



state energy for 4 He for both calculations as a func- 
tion of the oscillator frequency htu in small model spaces 
(e max =l, 2 and 3). The MCSM results obtained with 
32~50 MCSM dimensions are in reasonable agreement 
with the results from a full diagonalization in these model 
spaces. 

The treatment of spurious center-of-mass motion of 
4 He is illustrated in Fig. |3l Figure [3] shows the depen- 
dence of the expectation value of H c m , and the ground- 
state energy (inset) obtained with Lawson's prescription 
parameter /3 c . m . in the e max =l, 2 and 3 model spaces. 
The expectation value of H c . m . decreases rapidly and 
reaches a converged small value. In this way, the spu- 
rious center-of-mass motion can be suppressed to a large 
extent by choosing a suitable (3 c .m. value. 



IV. LOW LYING SPECTRA OF 10 Be AND 12 Be 

The 10 Be nucleus is a good candidate for testing ab ini- 
tio calculations employing realistic nuclear interactions, 
as there are adequate experimental data both in the 
ground state and in the excited states, e.g., excitation 
energies of two J 7r =2 + , T=l states and the B(E2) value 
of those states to the ground state. The AMD calcula- 
tions of Be isotopes [4l( , the GFMC approach @, Q and 
the NCSM [ll|, [42[ have been used to investigate p-shell 
nuclei like 10 Be and to reproduce features such as bind- 
ing energies and excitation spectra. This work is a new 
attempt to investigate these states by applying the no- 
core MCSM with realistic nuclear interactions. In this 
section, we present MCSM results for 10 Be and 12 Be. 

We discuss MCSM calculations using Vucom(N 3 LO) 
potential in an e max =3 model space. For the beryllium 
isotopes, the ground-state energies exhibit a minimum for 
oscillator frequencies huj around 16.0 MeV in the conven- 
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FIG. 4. (Color online) (a) Excitation energies of the 2f state 
(green diamonds), the l\ state (blue triangles) and the 
state (red squares) as functions of the MCSM dimension in 
the e ma x = 3 space, (b) Energies of the states shown in (a) 
as well as the Of (ground) state (black circles). f3 c .m. ■ fuv/A 
= 10 MeV is used to remove spurious components. 



tional shell model calculation [43j. We use bare charges, 
hence the electric quadrupolc moment is equal to the 
proton quadrupole moment. 

For a more precise investigation, we have to remove 
spurious components with respect to the center-of-mass 
motion, if they are mixed in calculated eigenfunctions. 
As discussed earlier, we use Lawson's prescription with a 
suitably chosen (3 c . m . in eq. (|5|). We shall use f3 c . m . -fuj/A 
= 10 MeV hereafter, unless otherwise specified. The 
same value has been taken in many MCSM calculations 
(not of ab initio type), e.g., [29| . 

The convergence of low-lying excitation energies as a 
function of the MCSM dimension has to be examined, 
as our goal is to investigate the excitation spectra. Fig- 
ure S] (a) shows the excitation energies of the 2^ , 2 \ 
and 0^ states as functions of the MCSM dimension for 
the e max =3 model space. If we evaluate the energy dif- 
ference e between results corresponding to the last two 
consecutive MCSM dimensions, we obtain e=18 keV for 
2+, 15 keV for 2+ and 37 keV for 0+ . The relative accu- 
racy of these excitation energies is ~ 0.3% for 2 + states 
and ~ 0.7% for 0+ state. In the MCSM calculation, the 
diagonalization is performed in a subspace comprised of 
25 to 50 optimally generated basis states. The size (di- 
mension) of this subspace is quite small compared to that 
of the entire Hilbert space taken in the direct diagonal- 
ization in the conventional shell model. This advantage 
will be even more obvious for heavier nuclei by the fact 
that the full diagonalization in e max = 3 is hardly feasible 
with other calculational techniques available presently. 

Figure |4] (b) exhibits the energies of the 0^~ 2 and 2+ 2 
states as functions of the MCSM dimension. One sees 
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steady improvements of these energies, particularly for 
the dimension greater than 30. The energies appear to 
become converged to a rather good extent. Figure |4] (b) 
shows that the ground-state energy becomes about -35 
MeV for dimensions large enough. This is still far from 
the experimental value ~ -65 MeV. After the Coulomb 
correction ~ 5 MeV, the difference is ~ 35 MeV, which is 
in part due to the choice of the interaction where three- 
body forces are missing. The model space and the con- 
vergence also contribute to the discrepancy. Such prob- 
lems are important issues in present and future MCSM 
projects. On the other hand, the present calculation ap- 
pears to be rather reasonable for excitation energies as 
shown later, and we use the Vucom(N 3 LO) potential in 
this first attempt. 

We now discuss properties of the and 2+ 2 states of 
10 Be. Figure [5] shows energy levels of these states. Some 
other low-lying states of 10 Be are shown also, and will 
be discussed later. While the MCSM results are about 1 
MeV and 0.6 MeV higher than the experimental values 
for 2f and 2 J, respectively, the basic patterns and scale 
are reproduced well by the MCSM calculation. In partic- 
ular, the low-lying 2^ level is a characteristic indicator 
of triaxial deformation, as discussed later. 

We now investigate these excited states in terms of 
the quadrupolc moments, E2 transitions and occupation 
probabilities. The quadrupole moments of protons and 
neutrons for the 2+ and 2\ states of 10 Be are shown, 
respectively, in Fig. [SJ One finds that beyond MCSM 
dimension of 30, those quadrupole moments reach stable 
values. The nucleus 10 Be has a negative quadrupole mo- 
ment for the 2f state. In contrast, the 2^ state shows a 
positive quadrupole moment. These features are also pre- 
dicted in Ref . [l| . We note that the protons have stronger 
deformation than neutrons in both states of 10 Be, be- 
cause there are two valence protons and four valence 
neutrons in the p-shell in major configurations, and the 
former produce stronger deformation than the latter. 

The B(E2) values from the 2+ 2 states to the ground 
state of 10 Be are shown in Tabled in comparison to re- 
sults by NCSM [Till!], GFMC 0, and AMD gj. The 
present result is rather similar to the CDB2k NCSM re- 
sults among those shown in this table. 

Some B(E2) values are calculated also for the mir- 
ror nucleus, 10 C, in the isospin formalism, as shown 
in Table HU This table indicates that MCSM value of 
B(E2;2 ] I " — > 0^ s ) appears to be in rather good agree- 
ment with the corresponding experimental data [44l [45[ 
for both 10 Be and 10 C. This is of certain importance be- 
cause from the viewpoint of the liquid-drop model, B(E2) 
value is proportional to Z 2 , and thereby the value of 10 C 
is expected to be larger than the corresponding one of 
10 Be, by a factor of 6 2 /4 2 in a naive expectation. While 
we take only bare charge (e p = e and e„ = with e being 
the unit charge), we can still produce almost the same 
values of B(E2;2+ -> 0+J of 10 Be and 10 C. This is be- 
cause although there are two more protons in 10 C than 
in 10 Be, they do not necessarily increase quadrupole de- 



formation, partly due to the OP3/2 closed-shell formation. 

We note that the B(E2; 2+ -> 0+ J value for 10 C has 
been obtained by NCSM calculations as 5.702 e 2 fm 4 [U 
and 10 ±2 e 2 fm 4 A GFMC value has been reported 

as 15.3 (1.4) e 2 fm 4 0. The present value, 9.3 e 2 fm 4 , 
appears to be the closest to the observed value. 

The MCSM value of the spectroscopic quadrupole mo- 
ment of the 2+ state of 10 C is obtained also from Fig. [5] 
as 3.04 e fm 2 by exchanging proton and neutron. 

The nuclei 10 C and 10 Be belong to the same isospin 
multiplct of T=l. In the notation of Timmer [47| . which 
makes direct use of the isospin formalism, one may write 
the E2 strength as 

B{E2) = [(ep + e n )S + T z (e p - e n )V] 2 , (7) 

where the e p and e n are the effective charges being e p = e 
and e n = in the present work. The reduced isoscalar 
and isovector matrix elements S and V must either be 
determined from experiment or be calculated with the 
help of suitable model wave functions. In Ref. [48|, the 
B(E2) value is proportional to [3.2 + 0.1 x T z ] 2 . The 
B(E2) value of 10 C should then be smaller than that in 
10 Be, as T z = -1 for 10 C and T z = 1 for 10 Be. 

Assuming that the Of and 2+ states of 10 Bc belong to 
the same K = rotational band, the intrinsic quadrupole 
moment Qo can be evaluated from the B(E2;0 ] f — > 2^~) 
value and the spectroscopic quadrupole moment with the 
following relations \^ 



Qo = 



(J+l)(2/ + 3) 
3K 2 - 1(1 + 1) 

^ • B(E2) t 
5 



Q, 



"I 1/2 



(8) 
(9) 



where Q is spectroscopic quadrupole moment, K stands 
for the K quantum number, and / is the angular mo- 
mentum of a member of the rotational band. The intrin- 
sic quadrupole moment evaluated by the spectroscopic 
quadrupole moment is 20.5 efm 2 , which is consistent to 
the one (21.6 efm 2 ) extracted from the B(E2;2f -> Of) 
value. This similarity seems to suggest an axially sym- 
metric deformation in the yrast band. On the other hand, 
the B(E2;2 2 f ~ — > 2+) is sizable, which hints at a notable 
triaxial deformation of 10 Bc. If 10 Be has strictly axial 
deformation, this B(E2) value should be hindered, as the 
transition between 2^ state and 2+ state is forbidden by 
the selection rule of K quantum number. The triaxiality 
leads to breaking of the K selection rule. For instance, 
B(E2; 2+ -> 0^)= 0.32 e 2 fm 4 leads us to a triaxial defor- 
mation with 7 = 11.4° in the Davidov-Fillipov model [50j |. 
Thus, the present results arc of interest in view of nu- 
clear shapes, although it may be an open question as to 
whether the classical picture of shapes can make sense in 
such light nuclei. It has been discussed in 5l| that 10 Be 
is triaxially deformed in a molecular-orbit calculation. 

Table E] shows B(E2;2+ 2+) and B(E2;2+ -> 0+) 
of 10 C too. These are considerably larger than the cor- 
responding values of 10 Be. At the first glance, the triax- 
iality appears to be more developed in 10 C than in 10 Be. 
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FIG. 5. (Color online) Some low-lying spectra of 10 Be in the e max = 3 model space. Black bars indicate experimental levels. 
Red levels are theoretical results obtained with the suppression of spurious center-of-mass motion (/J c . m . ■ Tiw/A — 10 MeV). 
Blue levels are obtained without removing the spurious center-of-mass motion. 



TABLE I. B(E2) values (e 2 fm 4 ) of 10 Be obtained by NCSM with CD-Bonn and CD-Bonn 2000 potentials [Tl|,H, GFMC with 
AV18 potential and AV18 plus different three-body forces 0, AMD HH, present MCSM, and experimental j4J,|i3| values. 





NCSM 


GFMC 


AMD 


MCSM 


Expt. 


Quantity 


CD-Bonn CDB2k 


AV18 AV18+IL2 AV18+IL7 








B(E2;2+ -+0+,.) 


6.58 9.8(4) 


10.5(4) 8.1(3) 8.8(4) 


9.46 


9.29 


9.2(3) 


B(E2;2+ ->0+,) 


0.13 0.2(2) 


3.4(2) 3.3(2) 1.8(1) 




0.32 


0.11(2) 



But one has to be careful, as this property holds for the 
proton part. As the proton part and the neutron part are 
exchanged between the mirror nuclei 10 Be and 10 C, it can 
be stated that the part consisting of four protons in 10 Bc 
tends to be deformed rather strongly in a prolate shape 
and the rest part (six neutrons) tends to be deformed 
in a triaxial shape, and the situation is just reversed in 
10 C. The deformation here may be static, dynamic, in 
between, or even molecular [5l[. Such a difference be- 
tween proton and neutron sectors is quite intriguing, and 
experimental investigations on these theoretical findings 
are of much interest. 

Figure [5] shows the 2+ 2 levels calculated in two ways (a) 
with center-of-mass motion suppression (default setup in 
this paper as stated already) and (b) without it for the 
sake of comparison. It is found that the suppression of 
the spurious center-of-mass motion keeps the excitation 
energies of the 2+ and 2\ states almost unchanged. It 
seems that the suppression of the cm. motion is not so 
relevant to these states. 

On the other hand, the cm. motion suppression is 
essential to the 0^", lj~ and 2\ states, as their excitation 
energies are raised by several ~ 10 MeV. 

It is likely that the lj~ state is sensitive to spurious 
center-of-mass contamination. We calculate the excita- 



tion energies of 1J~ state of 10 Be with j3 c , m .tiui / A = 0, 10 
and 20 MeV. The levels are shown in Fig. [7] When 
changing from f3 c . m huj/A = to 10 MeV, the excita- 
tion energy is increased by about 4 MeV. The result with 
Pc.m.ftw/ A = 20 MeV is less than 1 MeV higher, present- 
ing a hint at convergence, although it is slower than the 
2^ state. 

The Oj and 2^ levels in Fig. [5] arc lying quite high 
compared to the experiment. These states are expected 
to be intruder states with a large amount of 2p2h and 
higher excitations from the p-shell. The e ma x=3 space 
is considered to contain the major configurations of such 
intruder states, but a somewhat larger space- is needed to 
stabilize those configurations by coupling to even higher 
configurations. 

According to the adequacy of the e max =3 space, the 
states being discussed are divided into two groups : (i) 
0+ and 2+ 2 states, (ii) 0^ and 2jj~ states. For the present 
interaction and nuclei, the e max =3 space appears to be 
adequate for (i), whereas is still too small for (ii). The 
calculated 2j~ and 3j~ levels with f3 c . m hui/A = 10 MeV 
are higher than experimental ones in Fig. [5l although the 
deviations are smaller than those for 0j and 2^" states. 
It is likely that the properties of 2^ and 3^ states are be- 
tween the groups (i) and (ii), and should be investigated 
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TABLE II. B(E2;2+ -> 0+ s .), B(E2; 2+ -> 0+ s .) and B(E2; 2+ -> 2+) values (e 2 fm 4 ) of 10 Be and B(E2; 2+ -> 0+ s .) of the 
mirror nucleus 10 C obtained by the MCSM and the experimental data [il. |45|. 





iu Be 






B(E2; 2+ -> 0+ s . 


B(E2; 2+ -S- 0+ s . 


B(E2;2+ - 


,2+) 


B(E2; 2+ -S- 0+ s 


) B(E2; 2+ -> 0+ 3 . 


B(E2;2+ ^2+) 


Exp. 
MCSM 


9.2(3) 
9.29 


0.11(2) 
0.32 


3.28 




8.8(3) 
9.30 


2.15 


12.81 
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FIG. 6. (Color online) Spectroscopic quadrupole moments for 
protons (black circles) and neutrons (red triangles) as func- 
tions of the MCSM dimension. Upper (lower) panel is for the 
2+ (2+ ) state of 10 Be. 
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FIG. 7. (Color online) The 1J~ excitation energy of 10 Be 
with different /3 c . m . values compared to the experimental data 
(black bars). Blue, red and green bars indicate, respectively, 
results with p c . m .hu/A= 0, 10 and 20 MeV. 



also with wider model space. 

Figure [5] indicates occupation numbers of single- 
particle orbits for the ground and some low-lying states 
of 10 Be. Left side is for the results without the cm. mo- 
tion suppression, whereas right side is for those with it 
(default setup of this work). We begin with the group 



(i), for which left and right sides do not show much dif- 
ference. This is consistent with almost unchanged level 
energies of 2+ 2 between the two corresponding calcula- 
tions, as depicted in Fig.[SJ Figure [5] shows that protons 
and neutrons are mainly in the 0s and Op orbits. The oc- 
cupation number of the Osi/2 orbit is about 1.8 for both 
protons and neutrons. This value is remarkably constant 
for the three states in the group (i), and changes very 
little between proton and neutron. This means that the 
0si/2 orbit is occupied by ~90 % probability for both 
proton and neutron. 

If the Qs 1 / 2 orbit is fully occupied, the 4 He core is ide- 
ally formed. The present result suggests that the prob- 
ability of the 4 He-core formation is about (0.9) 4 ~ 2/3. 
The breaking of the 4 He core is nothing but the polar- 
ization of the core, which yields effective charges in the 
shell model with a core. In the present calculation, this 
effect is explicitly treated, producing the right amount of 
B(E2) values as discussed above. The UCOM transfor- 
mations act on short-range part of relative-motion wave 
functions. Electromagnetic operator at long wave length 
limit is then expected to be unaffected to a large extent. 
Thus, we use bare charges and E2 operator, for simplic- 
ity. 

The occupation number of the sd shell turned out to 
be about 0.7 for protons and neutrons combined, which 
corresponds approximately to the number of nucleons ex- 
cited from the 4 He core. More precise studies on the 
process of effective charges will be of much interest. 

For the group (ii) (0^ and 2jj~ states), the occupation 
numbers do change substantially between left and right 
sides of Fig. [5] This is consistent with the changes of 
their excitation energies shown in Fig. [5j 

The MCSM results for 10 Be low-lying spectra can be 
compared with those of NCSM in Ref. [nj|. Although the 
NCSM calculation does not use the same potential, the 
MCSM calculation shows similar results with the NCSM 
calculation for the 2^~ and l\ state. More results of the 
NCSM are listed and discussed in Ref. [ll[ for further 
comparison. The NCSM approach may have some dif- 
ficulty for similar calculations because the full sd shell 
configurations cannot be included at 8hu> truncation, for 
instance. In the present MCSM, the sd configurations 
are fully included. 

As we know from the experimental data, the 2f state 
of 12 Be has a lower excitation energy than the 2^~ state 
of 10 Bc. This is interpreted as a phenomenon related 
to the evolution of the magic number in exotic nuclei. 
However, the present MCSM results in an e max =3 model 
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FIG. 8. (Color online) The single-particle-orbit occupation numbers of the 0^ state (hollow columns), the state (full red 
columns), the 0^" state (full blue columns), the state (full green columns) and the 2jj~ state (full black columns) for 10 Be 
protons (upper panel) and neutrons (lower panel). The wave functions are calculated with f3 c . m .hui/A = (a) MeV and (b) 
10 MeV. 



space do not show this feature as can be seen in Fig. HO 
The suppression of spurious cm. motion may make the 
discrepancy to experiment larger. We definitely need a 
larger model space, and it is not tractable presently. 



V. SUMMARY 

For the first time, we have applied the no-core MCSM 
with realistic UCOM-transformed interactions to the in- 
vestigation of structure of 10 Be and 12 Be. 

We calculate some low-lying states of 10 Be and 12 Bc in 
an e max =3 model space. The results for the 2f and 2~£ 
states of 10 Be show a reasonable agreement with exper- 
imental data. We have kept particular attention to the 
spurious center-of-mass motion, and have suppressed by 
Lawson's method throughout this work. The MCSM re- 
sults show negative (positive) quadrupole moment for the 
2^ (2^~) state for 10 Be. We have analyzed the sensitivity 



of physical observables to center-of-mass contaminations. 
Some states, e.g., the 2\ 2 , are stable against variations 

of Pc.m.i but others, e.g., the lj~ state, are sensitive to 
fic.m.. In fact, the suppression of spurious cm. motion 
moves the lj" level by about 4 MeV to a better agree- 
ment to experiment. The obtained B(E2) values are 9.29 
c 2 fm 4 for B(E2; 2+ ->• 0+J for 10 Bc and 9.30 c 2 fm 4 for 
its mirror nucleus 10 C, which are close to the experimen- 
tal data with proton bare charge. Intrinsic quadrupole 
moments of the 2^ of 10 Be, evaluated from the spectro- 
scopic quadrupole moment and B(E2) value are similar to 
each other, suggesting an axially symmetric deformation. 
However, calculated B(E2; 2\ — > 2+) value is sizable, be- 
ing consistent with a modest triaxial deformation. The 
triaxial deformation is predicted to be more developed in 
10 C, providing intriguing issues on the mirror nuclei 10 Be 
and 10 C to be further investigated. 

The MCSM calculation presented here were performed 
in model spaces up to e max =3. Additional results, e.g. 
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FIG. 9. (Color online) The low-lying excitation energy levels 
of 12 Be calculated by the MCSM in e max =3 with the UCOM 
potential. The black bars indicate experimental data [52l |. 
while the red bars the MCSM results with /3 c . m .fk>j /A — 
MeV. 



for 10 Bc in e v 



=2, can be found in Ref. [42 



Al- 



though the ground-state energy of 10 Be is changed by 
about 10 MeV when going from e max =2 to e max =3, other 



observables turn out to be more stable already for the 
e max=3 model space. Thus, while the use of larger model 
spaces in the MCSM is certainly interesting and resultant 
changes will be evaluated in future, excitation energies 
and transition matrix elements look converged to a cer- 
tain extent. 
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